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Extending the Balitsky-Kovchegov (BK) equation independently to running coupling or to fluc- 
tuation effects due to pomeron loops is known to lead in both cases to qualitative changes of the 
traveling-wave asymptotic solutions. In this paper we study the extension of the forward BK equa- 
tion, including both running coupling and fluctuations effects, extending the method developed for 
the fixed coupling case [f[. We derive the exact asymptotic behavior in rapidity of the probabilistic 
distribution of the saturation scale. 
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In the leading-logarithmic (LL) approximation of high-energy QCD, the evolution with rapidity Y of the dipole- 
target scattering amplitude which appears in the dipole factorization of hadronic cross-sections, is given by the Balitsky 
hierarchy 0, or by the equivalent JIMWLK equation 0, 0, i, i, 0, i, S Q3, and in their mean-field approximation, 
CNJ ! by the Balitsky-Kovchegov (BK) equation 0, [lj], [12] ■ These equations describe the transition from the dilute regime 
£jry to the dense (saturated) regime of hadronic matter at weak coupling. The BK equation complements the linear 
Balitsky- Fadin-Kuraev-Lipatov (BFKL) [3, [3 (3 equation, whose solutions are exponentially growing, by involving 
a nonlinear damping term ensuring the taming of the growth, i.e. saturation. The Balitsky and JIMWLK equations 
include correlation effects in the dense regime. 

The BK e quat ion has been mapped [H, E3, Gil on t° the Fisher and Kolmogorov-Petrovsky-Piscounov (FKPP) 
equation (ljj, |20||, which admits universal traveling- wave solutions [2l], Q. This mapping has allowed to rederive 

i | some analytical results found in [H| , to find various additional results [H, [24], [25[ , and to give a simple picture of 

high-energy QCD. In particular geometric scaling [26|], observed in the HERA data, corresponds to a trave ling -wave 
Q-f structure of the dipole-target amplitude, and can be understood in term of solutions to the BK equation [27] . The 
q. geometric scaling is the scaling of the dipole-target amplitude 

' 

kx being the transverse momentum scale of the dipole, Q S (Y) the saturation scale, r g the scaling variable. 

A key feature of the FKPP equation is that its traveling-wave solutions are driven by the dilute region. Hence 
\Q . these solutions are extremely sensitive to discreteness or disorder effects in the dilute regime which can be added as 
a noise term to the FKPP equation. The exponential growth in the dilute regime tends indeed to enhance the effect 
, of the fluctuations. The obtained equation is called the stochastic FKPP equation (sFKPP). In particular, if one 
interprets the FKPP equation as a reaction-diffusion equation for a system of walkers (which we will call particles), 
the particle density can be arbitrary small forward the front of the traveling-wave solutions, even much smaller than 
the density corresponding to one particle in the whole system, while this very forward region drives the wave solution. 
| The FKPP equation is thus only a mean-field approximation of a reaction-diffusion process, neglecting the relevant 
fluctuations due to the particle discreteness . These fluctuations can be separated between two types, which have 
different consequences. First, most of them can be taken into account by cutting off the particle density when it is 
unphysically small. In [28|] , the modification of the traveling- waves solutions of the FKPP equation in the presence 
of this cut-off has been derived. The velocity of the wave is reduced w.r.t. the mean field case, and the shape of the 
front is modified in the dilute region. In the small cut-off limit, the solution is very slowly converging to the mean 
field solution, typically through inverse of logarithms of the cut-off. Second, in I], the effect of large fluctuations, 
not taken into account by the cut-off formalism, has been studied. Those fluctuations, which bring some particles far 
away forward the wave front, provide a random shift of the front. Hence, the position of the front becomes a random 
variable. 

This problem of sensitivity to the discreteness of the interacting particles has been raised in the QCD context in 
[H, [13, HH, [HJ for the simulations of onium scattering with unitarity corrections, and discussed in H,dl] where they 
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are associated with pomeron loops. This has led to conjecture that high-energy scattering processes in QCD may fall 
in the same universality class as a reaction-diffusion process (35j |. extending the mapping from BK to FKPP beyond 
the mean field approximation. The generalization of the Balitsk y a nd JIMWLK equations including pomeron loops 
have been investigated 1 through various formalisms [H, HE [HEjlEJlgl, S3, S3, S3, SE S3|. For example, the QCD 
evolution equation with pomeron loops has been formulated in [34 . 38] as a Langevin equation. If one neglects the 
non-locality of the dipole-dipole interactions and performs a coarse-graining approximation, the evolution equation 
with pomeron loops can be rewritten as a simpler Langevin equation 2 [34|, namely 



dyAf(L, Y) = a x(~d L )Af(L, Y) - Af 2 (L, Y) + ^Ka 2 s JV(L,Y) V (L, Y) 



(2) 



where N is the forward dipole-target scattering amplitude, L is related to the dipole momentum scale kr and some 
reference scale Qq by L = log(k^/Qo), K is a coarse-graining factor taking into account the approximation, and 
a = N c a s /ir. %(— 8l) is the LL BFKL kernel restricted to zero momentum transfer and zero conformal spin, with the 
eigenvalues x(l) = 2W(1) — ^(7) — V?(l— 7). 77 is a Gaussian white noise with the correlation 



( V (L, Y) r,(L',Y')) = ^5{L- L')8(Y - Y') , (3) 
a 

and the equation ([2|) should be understood with the Ito prescription 3 . Our goal will be to study this equation with 
running coupling instead of fixed. 

The main prediction of the fluctuations studies at fixed coupling is the existence of diffusive scaling at high rapidity. 
The fluctuations indeed seems to randomly shift the wave front. Hence, geometric scaling holds event by event but is 
no more valid in average, and is replaced by diffusive scaling, i.e. 



*y,g)-*(r,- 1 °'<y» ). (4, 

DY being the variance of the random saturation scale. Moreover, the averaging reduces the slope of the amplitude. 
The diffusive scaling has been conjectured in [35J, and derived in the small noise limit in [II], using the results of [l[ 
applied to the equation @. It has been also derived in the strong noise limit (49[. Numerical simulations indicates 
that this diffusive scaling holds in the case of a finite noise term [HO, HH . The diffusion coefficient D is a growing 
function of the strength of the noise. 

Almost all the above-mentioned studies have been performed at LL accuracy. But the extension of the saturation 
equations beyond LL level is one hot topic. A part of the Next to Leading Logarithm (NLL) corrections to the 
BK equation is already known [52l . l53l [54 . [HJ . In general, some of the NLL corrections contribute to the running 
coupling and some other to the kernel and to the nonlinear terms. If one knows them, one can resum higher order 
contributions to the running of the coupling, and obtain an equation containing the running coupling at the relevant 
scale and NLL kernel and nonlinear terms. The precise shape of the nonlinear terms is not relevant at high rapidity 
in the dilute regime, due to universal properties of the traveling waves solutions. It has also been shown [561 ] that the 
NLL corrections to the kernel are not relevant for the transition to the saturation regime at high enough rapidity in 
the running coupling case. By contrast, the replacement of the fixed coupling by the running coupling has a significant 
impact on the asymptotic solutions [It], [H, Hz|- I n particular, one has at large rapidity \ogQ 2 s {Y) oc VY, instead of 
logQ 2 (y) a y in the fixed coupling case. 

As the fluctuations and the running coupling seem to be both relevant for the saturation at high energy, they have a 
priori to be both taken into account 4 . The purpose of the present work is to study the solution of a model extending 
the BK equation, with both running coupling and fluctuations effects. 



1 The aim of these studies is not only to include fluctuations, but also to build a formalism for saturation with a symmetric treatment of 
projectile and target. 

2 The approximations leading from the full equation with pomeron loops to J2Jl are not well understood. Moreover, some doubt has been 
cast on the validity of the approximation for QCD [38, 45j. However, J5J is a useful toy model to investigate the fluctuation effects 
in QCD, and we will keep this approximation in our work. 

3 See for example [4r3l or [47ll for a review. 

4 A first attempt in that direction can be found in [581 . 



We explain in section [IT] how we implement the running coupling in @ to build our model. Then, we apply the 
method developed in [l[ to our model, implementing the cut-off in section fill Al and the large fluctuations in section 
IIIIBI Our final result is presented in section IIVI In concluding section [Vj we compare our result to the one in the 
fixed coupling case, and discuss the relevance of our model. 



II. FORMULATION OF THE MODEL 



Let us take into account the running coupling effects in the equation As mentioned in the introduction, the full 
NLL corrections to the BK equation are being calculated, but not yet available. However they are not all contributing 
to the asymptotic behavior of interest in the present paper [56| . Even the precise choice of the momentum scale for the 
running coupling will have no effect on our results 5 . In the following, we take the running coupling at the saturation 
scale everywhere in the equation namely 



a-^bpsiY) where b = and p s (Y) = \og(Q 2 s (Y)/A 2 QCD ) . C) 

With the assumption ([5|), the evolution equation ([2]) becomes 



b Ps (Y) d Y N{L, Y) = X (-d L W(L, Y) Af 2 (L, Y) + ^ "^^ff r,(L, Y) , (6) 
with (r?(L, Y) r)(L', Y')) = 4b p s {Y) 5(L — L')8(Y — Y') , 

where L is now defined by L = \og(k 2 -,/ Aq CD ). As p s (Y) is a priori unknown, one start with the assumption that 
(J6j) has indeed a traveling front solution, p s (Y) being the position of the front in L variable, and then check for 
consistency. Thus p s (Y) and Af(L,Y) can be both determined by ©. 

Let us start our analysis by considering the Ansatz p s (Y) = v(Y@) Y@ where v(t) has a finite and non-zero limit 
for t — > oo, t = Y@ being the effective time of the wave propagation. Writing the equation ^ in the t variable, 
one remarks that the only consistent value of (3 is 1/2 (as expected from the previous studies of BK equations with 
running coupling [l7l |23|). Then, the evolution equation writes 



b jp- 9tN (L,t) = x (~d L )N(L,t)-N 2 (L,t) + ^^^ u{L,t) , (7) 
with (v(L, t) v(L', t')) = S(L-L')8(t-t') , 
where t = y/Y, N(L, t) = Af(L, Y = t 2 ), u(L, t) = rj(L, t 2 ) / v / 2b v(t), and R = 2n7r 2 /N 2 b . 



III. DERIVATION OF THE SOLUTION 



A. Cut-off contribution 



The fluctuation term in is relevant only at very small density, namely for N(L, t) < R/v(t) t 2 . Hence, the shape 
of the wave front should not be modified at larger densities. However, in the standard BK equation case, the selection 
of the wave velocity happens forward the front, where the density is very small. One expects that the fluctuations 
modify the forward part of the front, and thus the velocity of the wave front. Following (28|, we replace the fluctuation 



5 A priori, the scales present in the problem are the size the parent and daughter dipoles. The first NLL parts calculated for the BK 
equation give hints to determine the running coupling scale [52t l53j|. except for the fluctuation term, which does not come from the 
typical fan diagram structure of the BK equation. However, only a rough estimate of the running coupling scale is needed, because the 
result of the present paper is quite robust with respect to a change of the running coupling scale, especially for the o? s coupling present 
in the fluctuation term. In order to check the independence of our result from this choice, the following analysis is done in the appendix 
fX]with the choice of the parent dipole momentum, i.e. the Fourier conjugate of the parent dipole size, as the running coupling scale. 




term in ([7|) by a cut-off, which corresponds to a deterministic approximation for the bulk of the fluctuations (see Fig. 
1). When the density is small, a fluctuation may indeed set the density back to zero where it stays since zero is a 
fixed point. 

We first study the linear part of the front interior region (see Fig. 1), defined by 

<JV(L,t)<l, (8) 



v(t) t 2 

where we can neglect both the nonlinear and fluctuation terms. In this region 



b -jp-d t N(L,t)~ X (-d L )N(L,t). (9) 

We now switch to the comoving frame of the front, replacing L by £ = L — v(t)t + £o and N(£,t) = N(L,t), £o 
being some parameter. Then 



[dtNfo t) - (v(t) + tv(t))dzm, *)] = x(-dt)Nfo t) . (10) 

The free parameter £o is fixed by imposing that £ = corresponds to the separation between the regions where the 
nonlinearity is relevant or not (see Fig. 1). Then, the linear part of the front interior is defined by < £ < Fo(t), 
Fo(t) being the position of the cut-off in the comoving frame. In [28|, the velocity has been found decreased by a term 
depending on the cut-off through logarithm. In our case, the cut-off is decreasing like t~ 2 because of the running of 
the coupling. Hence, we can expect that the velocity is decreasing through powers of logarithm of t. We thus drop 
the time derivative term in (fT0|) . which is decreasing like t~ x . We obtain the differential equation 



b -4^(v(t) + tv(t)) d 5 N(^ t) = X (-dt)N& *) (11) 



whose solutions are 



(12) 



where the 7t(i) are the solutions of the dispersion relation 



b - v (t)(v(t)+tv(t)) = ^ 1 



(13) 



As the fluctuation term is vanishing for t — > oo, j(t) and v(t) should converge to the critical values j c and v c of the 
BK equation with running coupling [171 ] , namely 



xilc) = 7c x'ilc) , and v, 



= / 2x(7c) 



(14) 



Following the method of [28J, we expand the dispersion relation (fT3|) around these critical values, and find 



- ^ tfft>(t) + 0((^(t)) 2 ) = (MO) 2 + Omit)) 3 ) , 

2 2~f c 



(15) 



where 5j(t) = j(t) — 7 C and Sv(t) = v c — (Sv > 0, since the cut-off is expected to reduce the velocity). We have 
supposed that Sv(t) cx (log£)~^ with some positive £, then t Sv(t) cx (logt) _< = _1 is a subleading term at very large t 
in the expansion (|15p. Hence 



S"/(t) = ±iir 



6v(t) 



where A 



7rV(7c) 
2b~/ c v c 



(16) 



We can rewrite at large t the solutions (fl2l) as 



(17) 



The shape of the front at small £, away from the cut-off, has to be the same as the one without fluctuations 



N(£, t) ~ A £ e" 7 ^ . 
One finds the constants C and D by matching (|17[) with (0, and writes 



(18) 



7T V Sv(t) 



a x / Mum \ „ , 

sin I W — e" 7 ^ for < £ < F (i) and <> 1 



(19) 



At the cut-off, one has 



N(F (t),t) = - 



TT y Sv(t) 



A 



«(t) 



(20) 



As Fo(t) and #u(t) are both unknown, we have also to ensure the continuity of the derivative of the front at Fo(t) 
in order to determine them both. The idea of the cut-off method is to prevent the front from having a too large tail 
in the dilute region. So let us impose at Fo(t) a stronger slope 7 > j c , with for example 7 — 7 C ~ 1. This gives a 
second equation 



IcA A 



5v{t) 



sin 7T 



^F oW ) + AcosL/M FoW 



-7c-Fo(i) = ^ K 

v(t) t 2 



(21) 



The equations (|20p and (|2"Tj) are similar to the ones in [28j, and gives the following result 

F 2 (<) Vc 2b lc v c F*(t) 



V W = v c- T^TTT = «c - ^7 rfTTT ( 22 ) 



AT(^,y) = sin ( I e~ 7 ^ for 0<f<F (t) and i large (23) 



7T 



*b(t) 



J e~ 7<=c for 0<£<F (i) and i large 



1 , / Aw(i)i 2 \ 2 , 1 , / Av c \ , ,, 

F i - - log F = - log* + - log c — + 0(Sv(t)) 

7c V(7-7cW 7c 7c V(7-7c)k/ 

These results have only a weak dependence in the free parameter 7 and in the coarse-graining factor n. In fact, as we 
have already truncated asymptotic series of logi in (|15p . only the first term is really significant. Thus, our result is 

2 

F (t) = —\ogt for the width of the front, (24) 

7c 

v(t) — v c 7c X (7c) velocity of the front. (25) 

8?w c log t 

N(£, Y) = 2Alogt S i n ( JHsL ) e -7e? for thc shape of the front; for < ^ < — logt and t large . (26) 

B. Contribution from large fluctuations 

As explained in [l|, the cut-off method allows to take into account the frequent but small fluctuations, which can 
eventually set N to zero. But the large fluctuations, which can give a non-zero value to N far away forward the front, 
are too rare to be modeled in a deterministic way. However they can have an important effect. The physical picture 
is the following (see Fig. 2): the traveling- wave is essentially given by the cut-off contribution, but sometimes, a 
fluctuation brings a few particles forward the tip of the front, setting there N of the same order as the cut-off. Then, 
this seed grows, forming a new forward-wave front, and also a backward-wave front. Finally, the old wave front meets 
the backward- wave front, and only the forward front created by the fluctuation remains. The long-time effect of a 
large fluctuation is thus to shift the wave front a bit forward. 

The first assumption in [l| is that such a fluctuation has a probability to happen during dt at a distance between 
5 and S + d5 forward the front 

V{8) d6 dt = de' 1 " 5 dS dt , (27) 

with C\ being some unknown normalization constant. This probability law depends on the shape of the front because 
the noise term in the sFKPP equation depends on the particle density, as in equation ([7]). We keep the assumption 
(|2"Tf in thc running coupling case 6 . 

Following [l|, one has to study the effect of one single fluctuation, i.e. the relaxation of the front, starting from 
an initial condition at t = to corresponding to the cut-off traveling-wave solution plus a seed at a distance S (see the 
first picture in Fig. 2). to is taken very large, so that Fo(to) is large. The effective width of the front interior F(t) 
is larger than the one for the cut-off solution Fo(t), starting at F(to) = fb(to) + 8. In the comoving frame (with the 
velocity l|25p). this relaxation is due to diffusion on the interval [0,F(to)]. Hence, the typical relaxation time is of 
order (^(io)) 2 . For the linear front interior up to the fluctuation, we take a similar Ansatz as in namely 

NIL, - Fit) O (^|^, J^) , (28) 

where v e (t) is an effective velocity of the front, such that 

v e (t ) = v(t ) and R(5) = lim [v e (t)t-v(t)t] ; (29) 



The assumption l|27| l has been verified numerically in the sFKPP case, and the physical arguments supporting it hold for the equation 
0. 




FIG. 2: Relaxation of the front after the creation of a seed by a large fluctuation, a) shows the initial condition with the 
cut-off front and the seed, b) the growth of the seed and c) the merging between the fronts, d) represents the final state of the 
relaxation, with a front similar the initial one (in dashed line), but translated by R(S). Note that the solution is drawn in the 
comoving frame with the velocity determined with the cut-off formalism. 



R{5) being the shift of the wave front induced by the fluctuation. As v e (t)t—v(t)t evolves only during the relaxation, 
and R(S) is of order 1, as 5, one has 

d t [v e (t)t-v(t)t] cx (F (t ))- 2 . (30) 
As 7 C ~ 0.6275, F(t) — Fo(t) starts of order 1, according to (|27|) . and relaxes to 0. Hence, one infers 

d t [F(t) - F (t)} cx (F (i )r 2 and then d t Fit) <x (F (io))~ 2 • (31) 
Inserting the Ansatz (|28|) in 

d t N(L, t) = X (-d L )N(L, t) (32) 

with the assumptions (|30|31|1 . one gets 

bv c d T G{y, t) - n 2 X "( lc ) G{y, r) = X "( 7c ) d 2 y G{y, r) , (33) 

where 

L-v e (t)t + £ £-(v e (t)t-v(t)t) 

y = 



F(t) F(t) 
t-tp 



(34) 



We take as boundary conditions 



G(0,r) = G(1,t) = 0, (35) 

because the solution has to match (|18|) at small y in order to take the nonlinear term into account, and y = 1 
corresponds to the position of the cut-off for the wave front generated by the fluctuation. The value of the cut-off 
gives a boundary condition for G too small to be relevant within the accuracy of our calculations. 

The equation (J33|), with (|33|) . is the same as for the FKPP case and for the BK equation with fixed coupling, after 
a rescaling of r. Hence, as in these cases, the first Fourier mode of G is constant, and the others decay exponentially. 
Calculating this first Fourier mode for the chosen initial condition, and matching it with the cut-off solution (|^|) 
translated by R(S), one finds 



fl(i) = 4 log ( 1 + w)' (36) 

where Ci is some unknown constant of order 1 . 

We turn now to the problem of the long time behavior of the solution, including several large fluctuations. (|36|) 
tells us that a fluctuation is relevant if e 1 " 8 is at least of order (F (t )) 3 - Hence, using ([77]) . one can conclude that 
the typical time between two relevant fluctuations is of order (Fo(i)) 3 . For t large enough, this time is much larger 
than the relaxation time of a large fluctuation, which is of order (Fo(t)) 2 . Hence we can admit that the wave front 
has enough time to completely relax between two large fluctuations. 

The calculation of the time derivatives of the cumulants of the position of the front done in |l| holds, thanks to this 
property, to ([T7|) . and to ([3^|) . It gives 



_ 7r 2 X " ( 7c ) 31ogF (t) 

2b lc v c F%(t) + 7 2 (F (t)r 



d t (p s (t 2 )) = v c - 



dt { P :(t 2 )) c = nK{n) for n>2, (37) 

7? +1 (F (t)) 3 



where the (.) means, in mathematical terms, the average w.r.t. the noise realizations, and is interpreted in QCD 
context as the average w.r.t. the target's color fields configurations. The index c denotes cumulants instead of 
moments of random variables. 

IV. FINAL RESULTS 

In the result (|37|) only the constant C\Ci remain unknown. We use the same trick as in [l[ in order to guess it. 
As explained in [l| , the expression of the velocity in the presence of the cut-off is still valid if one adds large 
fluctuations, but with a larger effective front width 

F eff (t) = F (t) + —logF (t). (38) 
Then, by matching with (|37[) the expression obtained, one can find C1C2. One gets 



dt( Ps (t 2 )) 



7r 2 x" (7c) , 37ry( 7c ) logi^t) 



2b lc v c F$(t) b^v c (F (t)) 3 
2u _ n\C{n)^x"{lc) 
7c" +1 ^c (F (t)Y 



Up to now, our results are quite similar to the ones at fixed coupling [l|. However, as the method used in sections 
IIII Al and IIII Bl is valid only if Fo(t) ^> 1, which is equivalent to logY" ^S> 7 C , we have only studied the asymptotic 
behavior of the saturation scale cumulants, above some very large rapidity Yq. Hence, integrating above Yq, 



( Ps (Y)) = ( Ps (Y )) + v c (VY- y/¥ a ) 



7c7r 2 x"(7c) 



Vy 



dt 



vTo log * 



37cT 2 x"(7c) /-^ df log log * 
log 3 i 



(Ps(Y))c = ( P i :(Yo))c + 



8bv c 
n!C(nW'( 7c ) 

,n-2 



To 



for n > 2 . 



(40) 



Two of the integrals in (|40[) can be rewritten using the Logarithmic Integral Li [59[ . Then our final result is 



(p s (Y)) = (p s (Y )}+v c ^VY-^Y^ 



7c7r 2 X // (7c) 
8bv c 



Li 



y 



logVF 



Li V^o 



logvyj 



37c7r 2 x"(7c) log log t 



{p n s {Y))c = {p n s {Y )) c + 



8bv c 



To 



log 3 * 



i6 7 rv 



Li 



y 



fogv/y 



(logVF)' 



- Li JY, 



lY* 



fx* 



\o%V% (log-s/yo)' 



for 7i > 2 . 



(41) 



We will consider also the following approximated form of our result. Assuming that the logs are not varying much 
and can be factorized out of the integrals in (j4H|) , one obtains 



( Ps (Y)) ~ ( Ps (Y )) + v c (W - y/Y?} 



TcTrVyTc) VY - V% 
8bv c log 2 VF 



37c^ 2 x"(7c) / at ^ log log ^ 



\/y- Vy c 

■ibv c 

n\ gn)^x"(lc) Vy~V% 

8^~ 2 bv c log 3 Vy 



log 3 VF 

:/>;(V), - .,^n;,;, - W " V , 3 for n>2. (42) 



On Fig. 3, the increase of the variance of the logarithmic saturation scale is plotted for the result (|4"Tj) and for the 
approximated result (|42|) as a function of the effective time y/Y/b associated with the traveling-wave dynamics at 
running coupling. Precisely, the plotted function is 

Var(Y) - Var(Y ) = (p 2 s (Y)) c - (p 2 s (Y )) c . (43) 

We have also plotted on Fig. 3 the result for another model, mixing fixed and running coupling, in which one takes in 
a at the scale Q S (Y) and a 2 for the fluctuation term at the scale Q s (Yq). Hence, it features the running coupling 
mean-field dynamics, but fixed coupling cut-off and fluctuations. It shows a linear rise of the variance of the front as a 
function of the effective time, as expected by analogy with the full fixed coupling case. Running coupling in the noise 
term strongly reduces the rise of the variance of the front. The truncation (|4"2"f of our result leads to significantly 
underestimate that variance, and thus should not be performed. 

The threshold Yo has to be adjusted, choosing a compromise between the precision needed and the phenomenological 
relevance of the rapidity range. We expect that the result (T4"Tj) catches at least the qualitative features of the full 
solution for Yq ~ 12 (which corresponds to consider that logyj ~ 2.5 is large enough compared to 7c ~ 0.63), and 
maybe down to Yq ~ 10. Hence, our result might be relevant for forward physics at the LHC. 

In addition, an implicit assumption has been done in our calculations. We have assumed that the wave front is 
formed down to the low density region where the noise term is relevant, so that it can drive the whole solution. The 
rapidity interval needed for the solution to form such a wave front is called Yf orm . Hence, we have an additional 
constraint on Yq which is Yq > Yf orm . In the present knowledge of QCD evolution equations, no reliable estimate of 
Yf orm is available in the running coupling case. It is yet difficult to estimate analytically (even roughly). It requires 
a much better understanding of the non-asymptotic behavior of the solutions of the BK equation with running 




FIG. 3: Increase of the variance of the logarithmic saturation scale Var(Y) — Var(Yo) = (pt[Y)) c — (ps{Yo))c, as a function 
of y/Y/b — y/Yo/b. Our result (|41[) is the solid line. The dashed line is the approximation (|42|) . and the dashed-dotted line 
corresponds to the result for a model with running coupling mean-field dynamics plus fixed coupling noise effects. We have 
taken here Y = e 2 ' 5 ~ 12. 



coupling 7 . Moreover, no numerical simulation of the full BK equation with running coupling and pomeron loops 
effects has been successful. 

In our running coupling results (|41[) . the increase of the cumulants has the same n dependence as in the fixed 
coupling result of Hence, the calculation of [II] is also valid in the running coupling case when the values of the 
cumulants at Yq are small compared to the increase between Yq and Y we have calculated. It happens for Y > Yq if 
Yq ~ Yf orm or generically for large enough Y. We can conclude that at high rapidity, the dipole-target amplitude is 
expected to show diffusive scaling 

tf{Y,k*)=Af(T d = L -b>PA , (44) 



with (p s (Y)) and {p 2 s (Y)) c given by (HJ|. 



7 A rough estimate relying on the asymptotic behavior of the solutions has been attempted in ref. [60], but it does not work in practice. 



V. CONCLUSION AND DISCUSSION 



In this paper we have studied a model extending the BK equation, and expected to include both the main features 
of fluctuations related to pomeron loops and of running coupling. Using the method developed in we have been 
able to calculate the high energy asymptotic behavior of its solution, and in particular the asymptotic behavior of 
the cumulants of the random saturation scale (|4ip . The effective time scale corresponding to the diffusive processes 
contributing to the front formation or relaxation is very short compared to the effective time scale related to the 
decrease of the coupling. That is why the method of [l|, built for fixed coupling, works also in the running coupling 
case. 

Let us compare the results obtained in the running coupling and in the fixed coupling case. Our result (|39[) is 
rather similar to the one of [l[ in the fixed coupling case, but some features are modified. The first difference is 
that the effective time is aY in the fixed coupling case and y/Y/b in the running coupling case. That difference is 
already present in the mean field approximation. The second difference is the logarithmic denominators in the running 
coupling case, which lower the rise of all the cumulants. 

The most subtle difference is in the validity range of the formula (p?9")) . In our case, it is valid only at high rapidity, 
whereas in the fixed coupling case, the corresponding formula is valid for any rapidity, provided that the cut-off front 
is formed, i.e. above Yf orm . Hence, in the fixed coupling case, one determines the cumulants completely. This result 
has been used in [II] to demonstrate diffusive scaling (j4|). It comes from a Gaussian approximation, which is valid 
because the magnitude of the cumulants is not growing too fast with the index of the cumulants. 

On the contrary, in the running coupling case, we have been able to calculate only the asymptotic behavior of 
the cumulants. At smaller rapidity, the front is forming, and this evolution is deterministic. This first regime is 
characterized by geometric scaling. Then, at Y — Yf orrn , when the front is formed up to the region of small N where 
the fluctuation term becomes relevant, the cumulants start to rise almost linearly, which breaks geometric scaling. 
Then, as the coupling at the saturation scale decrease, the relevant fluctuations becomes rarer. The rate of growth 
of the cumulants thus decreases and the evolution of the probability law of the saturation scale slows down, without 
stopping. Our result (|41j) concerns only this last regime. A priori, a large contribution to the cumulants may come 
from an intermediate region in rapidity, between Yf orm and Yq. It might delay the onset of diffusive scaling and lead 
to a non-scaling regime between the geometric and the diffusive scaling regimes. Hence, our result is not sufficient to 
predict if diffusive scaling holds in the physical domain in the running coupling case. Nevertheless, for large enough 
rapidity compared to Yq the contributions calculated in the present paper will be large compared to the cumulants at 
Yq, and allow us to predict diffusive scaling (|44[) in the high energy limit. Moreover, numerical simulations indicate 
that the diffusive scaling at fixed coupling is valid not only in the small noise limit, but also at finite noise. One can 
thus assume that in the running coupling case, the value of the cumulants at Yq or below should also be compatible 
with diffusive scaling. Hence our result (|4"Tj) 

• proves the diffusive scaling (as formulated in (|44[) 1 in the high rapidity limit in the running coupling case; 

• provides a hint that the diffusive scaling should hold for smaller rapidity, perhaps down to Yf orm , in the running 
coupling case. 

It is important to note that diffusive scaling can occur only at rapidities larger than Yf orm , whereas geometric 
scaling is valid at smaller rapidities. Hence, a large Yf orm could explain the success of the geometric scaling to describe 
the HERA data. Thus, one infers an experimental lower bound Yf orm > 9. If Yf orm < 12, it should be possible to ob- 
serve diffusive scaling (|4"4")) at the LHC. If Yf orm is larger, geometric scaling will remain valid in the LHC rapidity range. 

Let us recall and discuss the main assumptions made in the building of our model. 

• The choice of the running coupling at the saturation scale can be a priori a matter of debate. It is indeed not 
completely consistent, as the saturation scale is supposed to be dynamically generated. For that reason, the 
choice of the running coupling at the saturation scale in the BK equation is the only one giving an equation which 
can be exactly mapped onto the BK equation with fixed coupling. Except the precise value of the parameters 
and the fact that the effective time is \/Y in the first case and Y in the second, those two equations give the 
same results. But other choices for the running coupling scale, involving the parent and daughters dipole sizes or 
momentum scales, give solutions with a slow convergence to the critical asymptotic front and more complicated 
subasymptotic dynamics. It is a major issue to understand the solution at intermediate rapidity, and hence 
for phenomenological studies. But it makes no difference for our result, valid only at very high rapidity, when 
the wave front has been already formed. In the appendix [XJ the whole analysis is redone with the choice of 
the parent dipole momentum as running coupling scale. The results are exactly the same. One can understand 
this fact in the following way. The relevant region governing the whole solution is the one mentioned as the 



linear part of the front interior. Indeed, the softer modes are cut-off by the saturation, and the harder by the 
fluctuations. The size Fo(t) of this region is growing, but only as logi^, whereas this region as a whole moves 
with the front as s/Y. Therefore, for two arbitrary logarithmic scales L and L' in that region at large rapidity, 
(L — L')/L is always suppressed as an inverse power of Y. The appendix El provides an explicit example of 
this mechanism. Thus we can conclude that any choice for the running coupling scale involving the parent and 
daughters dipoles sizes or momentum will give the same result. Note that, as mentioned in the introduction, 
the other NLL contributions cannot be relevant asymptotically (56|. Hence, our result (jUJ can be considered 
as safe from NLL and higher order effects. 

• Concerning the fluctuations, our approximations are obviously the ones giving the Langevin equation (T2|). In 
particular, the off-diagonal Gaussian noise [381 ] is replaced by a diagonal one, and the dipolc-dipole scattering is 
taken local. These approximations are supposed to provide, at least qualitatively, the same type of fluctuations 
effects as the full pomeron loops equations. However that is not sure, as neither analytical nor numerical 
solutions to the full pomeron loops problem are available yet. 

As an outlook, it would be interesting to obtain analytical results at smaller rapidity than in the present work. 
However, it would probably require a different framework than the method of [1[ used here. It would also be sensitive 
to more details of the model. It would be useful to have a better understanding of the front formation in the running 
coupling case, in order to have a theoretical prediction for Yf orm . The main open problem remains however the 
understanding of the solutions to the full pomeron loops equations. 

Note that our results might be relevant outside QCD, as a study of a reaction-diffusion process with space-dependant 
reaction and diffusion rates. A chemical system with a temperature gradient could show such features. 

While completing this work, results of the numerical simulation f&oj ] concerning the extension with running 
coupling of a 1+1 dimensional reaction-diffusion model [5l[ for high energy QCD have appeared. The main result of 
(60j is that fluctuations effects are very strongly suppressed in the running coupling case, up to a quite high rapidity. 
Indeed, Yf orm is found to be large in that toy model, of the order of a few hundreds, and thus beyond the reach of 
experiments. It remains to study if such a large Yf orm is specific to that model, or occurs also in QCD. The results 
of [60l ] and of the present study seem however to be compatible concerning the asymptotic behavior of the cumulants 
of the saturation scale. Note also that the model used in [601 ] is a one-dimensional model with similar assumptions 
concerning the fluctuation term as the one studied here. 
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APPENDIX A: RUNNING COUPLING AT THE PARENT DIPOLE MOMENTUM SCALE 



In this paper, the running coupling has been taken at the saturation scale for convenience. But our results are 
more general, as they does not depend on the precise scale of the running coupling. As an example, let us present the 
derivation when one chooses the running coupling scale to be the parent dipole momentum fcx. Inserting 

fi5 4 (A1) 

in the equation ([2]) , one gets, instead of j6| , 



b L dyM{L, Y) = X (-d L W(L, Y) - Af 2 (L, Y) + ^ ^fJ^P r,(L, Y) , (A2) 

with (t](L, Y) t](L', Y')) =46 L 5(L-L')8{Y-Y') . 
Rewriting (IA2j) with t = \/Y instead of Y and with v{L, t) = r](L, t 2 ) y / t/2bL, one gets 



|| d t N(L, t) = X (-9l)N(L, t) - N 2 (L, t) + sj ^JJ^ v{L, t) , (A3) 
with {v{L,t)v{L',t')) =S(L-L')S(t-t') , 



in analogy with equation ([7]). The differences between the equations ([7]) and (|A3|) are more explicit if we use the 
comoving coordinate £ = L — v(t)t to rewrite (|A3|I as 



v(t) + y d t N(L,t) = x(-d L )N(L,t)-N'(L,t) + 
Hence, the cut-off is now for 

N(L,t) 



R N(L, t) 



\ t 2 (y(t) 



u(L,t) 



(A4) 



t 2 (v(t) 



Fo(t) 



(A5) 



instead of n/t 2 v(t). In the linear part of the front interior region, we have the equation 

b - (v(t) + [dtNfo t) - (v(t) + ti>(t))dtN(£, *)] = X (-dt)N(£, t) , 



(A6) 



instead of (fTU)) . In the next step of the method one neglects power suppressed terms. Thus one recovers the equation 
(fTTjl . Hence, the following steps are the same, up to the equations (|20[) and (|2"Tj) in which the explicit value of the 
cut-off appears. These two equations are replaced by 




sin tt\ 



I6v(t) 



Fo(t) 




Fo(t) 



-7cF (t) _ 



t 2 (v(t) 

-7cF (t) _ _ 



Fo(t) 



JK 



t 2 



(v(t) 



t 



(A7) 
(A8) 



The Fo(t)/t terms in the denominators of (|A7[) and (|A8[) gives a contribution to Fo(t) only at the order 0(Fo(t)/t). 
Hence, the results of the cut-off formalism (|25M24p remain the same. 

Let us now consider the study of large fluctuations. The equation describing the relaxation after one large fluctuation 
is again (|A6j) instead of (|32|) . The £ /t term is power suppressed in the linear part of the front interior, thus one finds 
again the equation ([55)1 . As the boundary condition for G at the cut-off is not sensitive to the exact value of the 
cut-off, the following results are the same. 

Hence, the results of our asymptotic study are identical if one chooses to take the running coupling either at the 
saturation scale or at the parent dipole momentum fcy. From the comparison of these two cases, one can infer that 
our results should hold for any reasonable choice of the running coupling scale. 
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